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' Abstract 

We consider Vela Jr. as being the old Supernova Remnant (SNR) at the 
^ . beginning of the transition from adiabatic to radiative stage of evolution. 

According to our model, Vela Jr. is situated outside Vela SNR at the distance 
of ~ 600 pc and its age is 17 500 yr. We model the high energy fluxes from 
^ Vela Jr. and its broadband spectrum. We find our results compatible with 

■ experimental data in radio waves, X- and gamma-rays. Our hydrodynamical 

. model of Vela Jr. explains the observed TeV gamma-ray flux by hadronic 

mechanism. The proposed model does not contradict to the low density 
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environment of the SNR and does not need extreme fraction of the explosion 
J> ' energy to be transferred to Cosmic Rays. 

O . Key words: Supernova remnants: hydrodynamics, gamma-rays. X-rays; 

Individual: RX J0852.0-4622 (Vela Junior) 
PACS: 98.38.Mz, 95.30.Lz, 95.85.Pw, 95.85.Nv 



^ ' 1. Introduction 



It is widely discussed that Supernova Remnants (SNRs) are among the 
most promising galactic candidates for the Cosmic Rays (CRs) accelerators 
up to the energy of 10^^ - 10^^ eV (see [H, 0] and references therein). The 
direct evidence for the acceleration of leptonic component of CRs in SNRs 
is synchrotron emission in a broad domain of the electromagnetic spectrum 
ranging from radio to X-rays. In the same time we are missing the direct 
evidence for the hadronic acceleration in SNRs. This is important point be- 
cause the composition of CRs is mainly hadronic according to the observed 
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CR flux and predictions of acceleration theories. The very promising mech- 
anism for the proof of the hadronic acceleration in SNRs would be the decay 
of neutral pions originated in the collisions of relativistic protons with pro- 
tons at rest. The 7r°-decay would result in 7-ray flux from the remnant or an 
SNR- molecular cloud system 0, 0, Isj. 

Therefore, after the detection of 7-rays from SNRs (with available instru- 
ments) it is crucial to know the nature (leptonic or hadronic) of the mecha- 
nism responsible for the production of the high energy flux. Both mechanisms 
inverse Compton (IC) and vr^-decay can contribute to the spectrum domain 
that modern instruments operate. The lack of a viable model that self- 
consistently explains physical conditions in the SNRs and lets unambiguously 
determine the production mechanism is important problem of the modern 
astrophysics and CR science. Even in well known detected by the H.E.S.S. 
collaboration 7-ray SNRs (RX J1713. 7-3946 0, 0, 8 and RX J0852.0-4622 
j^, 10]) the mechanism of the origin of 7-rays is still questionable. A broad 
MeV to multi TeV study of the objects is required to favour or reject the 



hadronic 7-ray production mechanism [11 



In our previous works [12, Il3| we showed that it is possible to expect the 



high energy 7-ray flux from the old SNRs at the transition from adiabatic to 
radiative stage of evolution. The high 7-ray flux from the SNRs evolving in 
either uniform or non-uniform ISM appears when the dense shell forms during 
the transition stage. The dense shell plays the role of the target material for 
the high energy CRs. We also showed that during the dense shell formation 
the thermal X-ray flux decreases because the significant amount of the SNR 
gas cools down. The low thermal X-ray and high 7-ray fluxes predicted 
by our model to exist in the same one SNR are the cases of the recently 
detected 7-ray SNRs (RX J1713.7-3946, RX J0852.0-4622). This gave us a 
hint to apply our model for the explanation of their nature. 

In the current work we apply our hydrodynamical model 0, H H [3 
of the transition stage of SNR evolution to explain the observed broadband 
spectrum (radio, X- and 7-ray) of Vela Jr. Supernova Remnant. The transi- 
tion stage itself was first studied numerically in 16|] and the estimates of the 
different reference times connected with the transition stage were reviewed 
in 



17 



Originally, Vela Jr. (or RX J0852. 0-4622) was discovered by Aschenbach 



in 1998 [18| from the ROSAT All-Sky Survey data. It is located in the 
south-eastern corner of Vela SNR, so in soft X-rays it is masked by the 
emission from Vela SNR and is visible only above ~ 1 keV. The emission 
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from Vela Jr. is dominated by non-thermal component with photon index 
r ~ 2.6 [l9|. The remnant is visible in radio band [2^. The radio flux is 
weak {Siy = 30 — 50 Jy at 1 GHz) with faint, limb-brightened emission similar 
to the X-ray morphology j2l|. The study of Vela Jr is complicated due to 
the high background created by Vela SNR. The distance and the age of the 
remnant are still unknown. The recent detection by Cherenkov telescopes 



22I . |9| of the TeV 7-ray emission from Vela Jr. and the fact that it is one 
of the few SNRs with non-thermal dominated emission made it an object of 
intensive studies 23, 19, 24, lol . 



In what follows we briefly describe the hydrodynamics of the transition 
stage (Section [2]). Then we apply our hydrodynamical model to Vela Jr. 
(Section [2D and derive the hydrodynamical and CRs parameters for the spec- 
trum modelling (Section H]). The broadband spectral model and 7-ray surface 
brightness map of Vela Jr. are presented in Section [51 We make a short dis- 
cussion of the obtained results in Section [6] and conclude in Section [3 

2. Transition Stage Model 

2.1. Origin and Dynamics of the Thin Shell. 



The full description of the transition stage model is given in [ij, [15|| . Here 
we outline the most important features. 



Numerical simulations [16|, [25|| show that radiative losses become impor- 
tant at the time: 

ttr = 2.9xlO'Eti''nfs!i'yr, (l) 

where E^i is the explosion energy in 10^^ ergs, tiism is the ISM number 
density in cm~^. 

Radiative losses lead to rapid formation of the cold dense shell near the 
front. When the hot SNR gas joins the inner boundary of the shell it cools 
down leading to the growth of the shell. The shell grows as well because it 
sweeps up the ISM. The transition phase ends when the hot gas stops cooling 
effectively and no more replenishes the shell. We use the numerical results 



16| saying that cooling is important for the hot plasma within the outer 
flve percent {a = 0.05) of the SNR radius at the beginning of the transition 
phase: Ar = aRtr- Parameter a is the only free parameter of our model of 
the transition stage of SNR evolution. 

For the beginning of the transition stage we take the time tfr- It is the 
time when the flrst cold gas element of the shell appears at the shock front. 
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The element has temperature Tgh, pressure Psh, density psh and velocity Vgh- 
From the balance of external and internal pressure on the shell pisMVsh — 
Ptriytr — Vsh)'^, we derive the velocity of the shell for the adiabatic index 
7 = 5/3: ^ 

Vsh = -Dtr = const (2) 

where pisM is the ISM density, ptr and Dtr are the shock front plasma density 
and velocity for the time Ur, Vtr = 0.75Dtr is the plasma velocity just behind 
the shock front. We assume that during the transition stage a small pressure 
gradient inside the SNR results in conservation of the velocity of each plasma 
element unless and until it joints the shell: 

.(a,t) = |y/"''-) (3) 
^ ' ^ [Vsh if a^t) <a<Rtr ^ ^ 

where adt) is the Lagrangian coordinate of the gas element that cools at the 
time t. For the end of the transition phase adtsf) — 0.78. The duration of 
transition phase in our model is: 

. CtRfr , ,\ 

2.2. Hot Gas Parameters Inside the Shell 

For the time ttr < t < tsf the velocity of the gas element with the 
Lagrangian coordinate < a < ac(t) is given by Eq. ([3]), so for the Euler 
coordinate r{a,t) we have: 

^\ i ^((^^ttj) +v{a,ttr){t -ttr) ifO<a<ac(t) 
' \ Rsh if ac{t) <a<Rtr ^ ^ 

The density distribution p(a, t) we find from the continuity condition: 

r(a,ttr) V dr{a,ttr) 



p(a,t)=p(a,t,.)(-^J (6) 



the hot gas pressure is: 
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and the temperature: 

where Rg is the absolute gas constant, is the molar mass. Thus we give 
the full description of the hot gas inside the SNR. This gives us a possibility 
to model thermal X-ray emission of the remnant. 

2.3. Cold Shell Gas Parameters 

Starting from the time ttr the mass of the shell increases because the 
interior hot gas cools and joins the shell and the shell sweeps up the ISM. The 
temperature of the shell equals the ISM temperature: Tgh = Tjsm = 10^ K 
and its pressure Psh equals the dynamical pressure on the shell. See Section H] 
for detailed discussion of the shell parameters. 



3. The Hydrodynamical Model of Vela Jr. 

In the current section we apply our model of the transition stage to 
Vela Jr. The high energy flux from Vela Jr. is difficult to explain by vr^- 
decay mechanism because in order to obtain the sufficient amount of p-p 
collisions the number density of the ISM in the vicinity of Vela Jr. should be 
high. The high density of ISM contradicts to very weak thermal component 
of the observed X-ray flux. The weak thermal X-ray emission puts an upper 
limit on the ISM density. The other possibility to explain the observed 7-rays 
by 7r°-decays is to have a higher fraction of the explosion energy transferred 
to CRs during the acceleration history. However, in this case we need almost 
all the explosion energy to be transferred to CRs (see Table 1 in which 
is not a plausible solution of the problem. 

As was shown above (Section [2]), the advantage of our hydrodynamical 
model is that the high density cold shell forms at the front region of the 
remnant during the transition stage. The CRs contained in the cooled volume 
of the SNR start interacting with the dense target material of the shell, thus 
providing an increased 7r°-decay 7-ray flux. We use this advantage to build 
the model of Vela Jr. 

We assume that Vela Jr. is rather old SNR that has gone through the 
adiabatic stage of evolution and is at the transition stage. We assume it is at 
moderate distance and do not consider a close-by scenario for several reasons. 
Firstly, because the observed absorption is rather high (N^ ~ 6.2 x 10^^ (lo|). 
Secondly, in our model it is old SNR and unless the density of ISM is high (but 
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this contradicts the observations) the radius should be quite large. Finally, 
there are some indications that Vela Jr. is not the close-by and young SNR 
(see Discussion Section for details). 

We have made the calculations of hydrodynamical parameters of Vela Jr. 
for the grid of models of SNRs at the transition stage. We stopped at the 
model that does not contradict to the very low thermal X-ray flux that 
is invisible in Vela SNR foreground and that shows very minor features in 
the harder band of the X-ray spectrum. Besides, there should be present 
rather high 7-ray flux due to 7r°-decays. Fitting the thermal X-rays from 
the model taking into account the absorption (we use Nh = 6.5 x 10^^) to 
the given observational upper limits we were able to reconstruct the energy 
of explosion, ISM number density and the radius of the SNR. The latter 
provided us with distance to the remnant. 

The model of Vela Jr. that we consider has the following parameters. 
The number density of the ISM at the place of explosion was n = 1.5 cm~^. 
The explosion energy was Esn = 0.2 x 10^^ ergs that is in range of the 
classical explosion energies of Supernova. The radius of the considered model 
is 10.2 pc that corresponds to the distance to Vela Jr. of about 600 pc given 
the angular size of 2°. Given that the transition stage starts at the time 
ttr = 16 000 yr, to explain the observed non-thermal fluxes (see Section [5]) 
the transition stage should have started 1 500 yr ago, so the considered age 
of Vela Jr. is 17 500 yr. 

4. Cosmic Rays in Vela Jr. 

Before the broadband spectrum modelling it is important to know the 
properties (i.e. the energy of CRs, their distribution inside the remnant and 
their spectrum) of relativistic particles that take part in different radiation 
processes. We must note here that we do not consider the physics of CR 
acceleration and its history in the remnant. We assume that the acceleration 
stopped at the end of Sedov stage when the SNR entered the transition stage. 
By this time, a fraction u of the explosion energy Esn was transferred to CRs. 

4.I. Cosmic Rays in the Hot Gas 

For the end of Sedov stage in the SNR there is a population of relativistic 
protons with the spectrum Jp{Ep) containing the energy uEsn'- 
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where Wcr is the energy deposit of CRs. The spectrum of particles is as- 
sumed to be a power-law with exponential cut-off: 



J,iE,) = A ■ E;^exp (10) 

where A is the normalizing constant, a is the power-law index and Ep^^ut is 
the cut-off energy of the proton component of CRs. 

For the primary electron component, we assume that it deposits a fraction 
Kep of the energy of the proton component and for the end of the adiabatic 
stage has the spectrum Je(-E'e) with cut-off energy Ee^cut- 

UE,) =K,p-A- E^'^exp ( (11) 

During the transition stage the electron spectrum evolves because of energy 
losses that are mainly due to synchrotron losses. Thus, to find the electron 
spectrum for the time t > tn- we solve the equation of the particle's energy 
evolution. We follow the procedure of how to find the spectrum of particles 
in the SNR for the given time described in [26| . 



In |12l . Il3| we investigated the case of uniform CR distribution within the 
SNR that is a good approximation for the end of the transition stage when 
CRs had enough time to diffuse into the inner layers of the remnant. Here in 
the case of Vela Jr. (the beginning of the transition stage) we adopt that the 
distribution of CRs repeats the distribution of the gas inside the remnant. 

4-2. Cosmic Rays in the Shell 

The shell formed during the transition stage is of special interest because 
according to our model it should provide the significant amount of vr^-decay 
7-rays. Based on the above assumption about CR distribution, the amount 
of CRs confined in the shell is given as: 

WcR.sh = ^TiRi [ u{r,ttrydr (12) 

Jr{acit),t) 

where u!{r,ttr) is the energy density distribution of CRs inside the SNR. 

The additional energy gain of relativistic particles can be considered due 
to adiabatic invariant conservation 27 . Adiabatic invariant is conserved 



when magnetic field changes spatially on the length scales much greater than 
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the particle gyroradius, Vg, and temporally on timescales much longer than 
the particle gyroperiod, Tg |28|, l29| . If we consider the CR proton of energy 



200 TeV gyrating in the magnetic field ~ 20 /xG, the particle would have 
Vg ~ 0.01 pc and Tg ~ 0.05 yr, implying that in order to apply the condition of 
adiabatic invariant conservation in cooling region of Vela Jr., the turbulence 
scales should be assumed to be much greater than provided values of Vg and 
Tg. The X-ray observations of Vela Jr. do not show any clumpy structure, but 
only smooth filaments located at the outer region of the shock [2j] . The lack 
of clumpy structure implies that there is no serious magnetohydrodynamic 



turbulence in the emitting plasma [30|, iMl]- The filaments are straight with 



length comparable to the radius of the SNR and widths ~ 50 arcsec (~ 
0.15 pc adopting our model distance) at shock downstream side 0]. If there 
were significant turbulent flows in Vela Jr. the filaments would be distorted. 
Finally, the fact that no time variability of the X-ray flux from Vela Jr. 
was reported is another argument against presence of any strong small-scale 
turbulence. Therefore, using condition of adiabatic invariant conservation, 
p\/B = const, where p± is the normal to the magnetic field B component of 
the particle's momentum, so the energy gain is ^JB^B2 = {psh/^PisuY^^ 
taking into account that B ~ p^^^. Therefore, the energy of CRs contained 
in the shell will be WcR.sh{psh/ '^PismY'^ ■ 

The spectrum of the CR protons will be almost the same as in the interior 
of the SNR, however with higher cut-off energy due to additional energy gain. 
No serious losses affect the proton component. 

The spectrum of the CR electrons is strongly affected by the losses, es- 
pecially synchrotron losses. Every new portion of the electrons coming from 
the hot interior of the SNR to the shell will effectively radiate in the high 
magnetic field of the shell [B ~ p^^^, so magnetic field in the shell is higher 
than in the SNR interior). Therefore, to calculate the spectrum of electrons 
for the time t > ttr in the shell we again solve the equation of the particle's 
energy evolution, but here we add the source of electrons that come from the 
hot SNR interior and join the shell. 

The energy of CRs contained in the shell is almost fully defined by the 
energy of CR protons. Having the energy of CRs we can determine the other 
parameters critical for the calculation of the broadband flux from the shell 
(i.e. number density of target particles, magnetic field). The shell pressure 
is given by the sum of CR, gas and magnetic pressures and equals to the 
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dynamical pressure on the shell: 



PCR + P,as + Pb = ^ + ^^^^ + ^ = PISmVI (13) 

o fi on 

where Ugh is the energy density of CRs in the shell. The dynamical pressure is 
determined from the hydrodynamical model of Vela Jr. We consider the case 
when there is an equihbrium between CR, gas and magnetic pressures. In the 
interior of the remnant thermal pressure dominates, but when the gas in the 
cooling region compresses and forms the shell the other pressure components 
become non negligible. According to ^3^, the density when magnetic pressure 
becomes important is Ucru = 0-7GVshn^^'^{BisM/ (J^G)~^ that in our case is 
~ 42 cm~^. Therefore, starting from this density in the shell the magnetic 
pressure is of order the gas pressure. If we considered only ISM CRs, then we 
would ignore CR pressure component because their ISM energy density is not 



enough to make significant pressure [32j. But we consider the compression of 



CRs contained in the SNR cooling region, where the energy density of CRs is 
high [33|. We assume that CR energy density would not exceed the magnetic 
energy density, because CRs are confined by the magnetic field and they 
would leak out of the shell until their energy density roughly equals magnetic 
one, therefore it is reasonable that all pressure components are roughly equal. 
From (fT3ll we obtain the number density Ush and the magnetic field Bgh of the 
shell as well as the fraction ly of the explosion energy transferred to CRs. If 
PcR is greater than other pressure components, i.e. u is greater than assumed 
in the current model, ngh and Bgh would be smaller. This would result in 
less efficient energy losses of the electrons and hence harder X-ray spectrum 
of the shell. If Pcr is smaller than other pressure components, the X-ray 
spectrum of the shell would be softer. In the 7-ray domain the 7r°-decay flux 
would not really change as decrease(increase) of the target protons would be 
compensated by the increase(decrease) of CR energy density. 



5. Spectral Modelling of Vela Jr. 

The observed broadband spectrum of Vela Jr. can be well explained by 
the combination of synchrotron radiation (radio, hard X-rays), weak soft 
thermal X-ray component, 7r°-decay 7-rays and some contamination from IC 
7-ray component. Each type of the emission comes from both the interior 
of the SNR and from the shell, so the final broadband spectrum consists of 
many components. To calculate the 7-rays and the secondary products of 
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p-p interactions we use the parametrizations of l3j] for the proton energies 
> 100 GeV and (5-function approximation from [35| modified in 3J] for the 
lower energies. 

5.1. Soft Thermal X-rays 

As our method is based on the hydrodynamics and we used thermal X- 
rays to constrain our model, we started from the thermal X-ray component 
that naturally displays the thermodynamic characteristics (density, temper- 
ature and pressure) of the gas flow inside the remnant. There is no thermal 
X-ray emission from the shell because it is cold. We calculate only contin- 
uum thermal X-ray emission as we are mostly interested in the energy range 
above 2 keV, where Vela SNR contamination is not significant. The lines 
contribution is known to be small in the energy domain above 2 keV. The 
thermal spectrum of Vela Jr. is shown at Fig. [H 

Let us explain two lines showing thermal spectra presented at Fig. [H The 
dash-dotted line represents the whole volume integrated absorbed spectrum 
of Vela Jr. model. However, we noticed that the analysis of X-ray data in 
Tol has been done using the regions observed by ASCA (Nl, N2, N7). 



The detected emission in the case of thermal component is defined by the 
integral along the line of sight within the volume projected to the observed 
region. Depending on the position of the observation region on the plane 
of the image, different SNR layers contribute to the observed emission. As 
known from the Sedov solutions the hottest layers are the innermost ones 
and the coldest are the outermost, so the thermal spectrum of the inner 
regions of the SNRs is harder than the spectrum of the outer regions. If we 
calculate the thermal spectrum of the ASCA regions covering the SNR shock 
by integrating from ~ 0.6-R to R, where R is the shock radius of the SNR, we 
obtain the spectrum presented by the dotted line at Fig. [H This spectrum 
is softer because the contributing layers are colder. It worth noting that this 
situation corresponds to the ASCA observations and the obtained spectrum 
(dotted line) is in good agreement with the spectrum presented in m. The 



upper limit of 1/30 of the X-ray flux for the thermal component in the energy 



range 2-10 keV put in [lOj is also well reproduced. 

The soft thermal X-ray flux (with estimated contribution of emission in 
lines on the level of ~ 5 times continuum emission) from Vela Jr. model 



here and below we mean thermal continuum component only 
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Figure 1: The absorbed X-ray spectrum of Vela Jr. in 0.5 - 10 keV range (solid black) as 
the sum of thermal " region" integrated (dotted green) and non-thermal (dashed red) com- 
ponents. The dash-dotted blue line represents the volume integrated thermal spectrum. 
The dark grey area shows ASCA data [3, 113] and the light grey area shows energy range 
contaminated by Vela SNR. 



obtained by integration over the whole SNR volume is consistent with non- 
detection of Vela Jr. at low energies on Vela SNR foreground. The flux from 



Vela SNR in range from 0.1 to 2.4 keV is 1.3 x 10 ^ erg/cm^s 36|] which is 
~ 18 times higher than Vela Jr. flux in the same range. 

5.2. Synchrotron emission 

The synchrotron radiation is responsible for the radio and X-ray flux from 
Vela Jr. There are two main components of the synchrotron emission: from 
the SNR interior and from the shell, both coming mainly from the primary 
electrons. In each of the components there is a small contribution from the 
secondary electrons born in 7r+/~ decays. The charged pions are born in 
collisions of relativistic protons with protons at rest in the same amount as 
neutral pions. The secondary electrons are taken into account by adding 
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the source term to the equation of the particle's energy evolution when it 
is solved to find the spectrum of electrons for the given time during the 
transition stage. The contribution from the secondary electrons is negligible 
for the interior, whereas for the shell it is noticeable, but still is very small. 
There are differences in the "interior" and the "shell" emission components. 
The former has the cut-off in softer range (< 2 keV) with a steep slope, the 
latter component's cut-off is in the harder band (> 2 keV) with a flatter 
slope. This is due to electrons constantly coming from the hot interior and 
radiating in the higher (compared to the interior) magnetic fleld of the shell. 
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Figure 2: The broad radio to X-ray unabsorbed non-thermal spectrum (soHd black) with 
contribution from the "interior" primary (dotted red), secondary (dash-dotted red) and 
the "shell" primary (long dash-dotted green), secondary (double dotted green) leptons. 
The radio data is taken from 21 1 , the dark grey area shows ASCA data [19|, |lO[ and the 
light grey area shows energy range contaminated by Vela SNR. 



In the current model we assume that the magnetic fleld in the interior of 
Vela Jr. is B = IS/xG and in the shell is Bg^ — Q7fiG that comes from the 
equilibrium of CR, gas and magnetic pressures. We consider that Vela Jr. is 
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1 500 yr in the transition stage. The cut-off energy of the electrons for the mo- 
ment tfr (the beginning of the transition stage) is assumed to be Ee^cut = 60 
Te V. To fit the X-ray flux to the observed value one needs the electron / proton 
ratio Kep — 0.016. Under the described conditions, for the simultaneous ex- 
planation of the radio and X-ray observations the power-law index of the 
electron component is taken to be a = 2.1. The X-ray synchrotron flux 
from Vela Jr. taking into account the absorption is presented at Fig. [H The 
X-ray flux in 2 - 10 keV band is Fx — 8.4 x 10~^^ erg/cm^s and is con- 



sistent with one obtained in |10l . Il9l | . The corresponding synchrotron radio 



flux is also consistent with experimental data (30 - 50 Jy) presented in [21 
The unabsorbed broad radio to X-ray range with contribution from different 
components ("interior", "shell") and the secondary electrons (positrons) is 
given at Fig. [2l The spectral index of the modelled X-ray radiation, F, in 2 
- 10 keV range is well reproduced. We derive F ~ 2.68. It is consistent with 



determined F = 2.79 ± 0.09 in Jlfl] and F = 2.6 ± 0.2 in [19 



5.3. The J -rays 

The detailed study of the 7-ray emission of Vela Jr. is presented in 
lo| . The observed flux in 1-10 TeV band is ~ 5.4 x 10^^^ erg/cm^s with 



photon index F ~ 2.24. Based on the fact that ISM number density cannot 
be high, the detected flux could be explained by vr^-decays only under the 
assumption that large fraction of the explosion energy was transferred to the 
energy of CRs. In our hydrodynamical model of the SNR evolution applied 
to Vela Jr. we use the advantage of the increased number density in the 
newly formed cold shell. The density in the shell increases by a factor ~ 37 
[ush = 55.2 cm~^), so we are able to explain the observed 7-ray flux from 
Vela Jr. with more conservative energy deposit of CRs that is ~ 4 % of the 
explosion energy. The parameters of the proton component are F = 2.15 and 
Ep,cut = 400 TeV. 

The 7-rays in our model are explained by vr^-decays. However, in the 
range 1-10 TeV we have IC contamination on the level of ~ 20% of the flux. 
The IC emission comes from the interior of the remnant where the magnetic 
field is not high enough for the IC emission to be suppressed. This is natural 
for the assumed normal conditions of the ISM in Vela Jr. vicinity. The high 
magnetic field in the whole Vela Jr. and its vicinity that would make IC 
contribution negligible seems quite unreasonable unless the SNR evolves in a 
very special region. However, in the shell formed during the transition stage, 
the high magnetic field is reasonable because of density increase {B ~ p"^^^)- 
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Figure 3: The total broadband spectrum of Vela Jr. (solid black) together with the exper- 
imental data (radio 21 1, X-ray 3, 7-ray IC contribution from the "interior" of 
the remnant is shown in dash-dotted green. 



It must be pointed out that without the transition stage with normal ISM 
density it is difficult to claim that the 7-ray emission comes from the 7r°- 
decays unless very high magnetic field inside the whole SNR and a large 
fraction of the explosion energy contained in CRs are assumed. 

The derived 7-ray spectrum perfectly fits the observed one of Vela Jr. 
The broadband spectrum of the remnant is presented at Fig. [31 

As was pointed out in Section the shell formed during the transition 
stage is thin. So, besides the explanation of the 7-ray fiux our model pro- 
vides a good explanation for the observed 7-ray morphology of Vela Jr. The 
surface brightness map of Vela Jr. model along with the H.E.S.S. image 
is presented at Fig. HI 
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Figure 4: The 7-ray H.E.S.S. data [9|(left) along with the model surface brightness map 
(right). The model was smoothed with the Gaussian of 0.1° to simulate the H.E.S.S. point 
spread function (PSF). The H.E.S.S. PSF is indicated by a circle. The contour lines are 
the ROSAT All-Sky Survey data above 1.3 keV. 



6. Discussion 

The breakthrough of modern observational techniques provided the com- 
munity with unprecedented results. The detection of 7-rays from long de- 
bated SNRs gives unique opportunity to find the clue to the solution of the 
high energy CRs problem. Nonetheless, the broadband data on the 7-ray 
SNRs appeared to be somewhat ambiguous in the sense of explanation of the 
observed 7-ray emission by the mechanism of the hadron origin. 

In this paper we aimed to give the explanation of the Vela Jr. broadband 
emission (radio, X-rays and 7-rays) applying the hydrodynamical model of 
the transition from adiabatic to radiative stage of the SNR evolution. Con- 
ditions that are formed during this stage satisfy the origin of hadronic 7-ray 
flux and thus can prove the existence of acceleration of CRs in SNRs. To this 
end we are comfortable with our result because it does not contradict to the 
experimental data and is self-consistent. A new age of 17 500 yr proposed 
here is the direct consequence of the hydrodynamical calculations for the 
normal ISM number density and assumed energy of explosion. The obtained 
radius of 10.2 pc and corresponding distance of 600 pc are well in range of 
the proposed earlier values (200 - 1 000 pc). Therefore, we find our result 



15 



for the age of Vela Jr. normal. The only thing that would contradict to 
this age is the Ti^^ line detection S^] which significance was questioned in 



38l |. Besides, there are strong implications that Vela Jr. is indeed rather old 
remnant. Among these implications is the measurement of the X-ray expan- 
sion of the north-western rim of Vela Jr. [39]. The expansion measurement 
showed that it is about five times lower than in the historical SNRs and is 
comparable to that of the Cygnus Loop which is known to be in the post 
Sedov stage of evolution [sP]. Another implication that Vela Jr. is not a 
young remnant is a possible association of the pulsar PSR J0855-4644 with 



Vela Jr. [40|. In case this pulsar is really associated with Vela Jr. the age 
of the remnant must exceed 3 000 yr or even more. Yet another object SAX 
J0852. 0-4615 could be associated with Vela Jr. Being located at the geomet- 
rical centre of Vela Jr. it is a neutron star candidate with a probable age of 



a few 10^ yr 41 



We would like to note that the fraction of SN energy transferred to CRs 
in our model of Vela Jr. is very near the conservative value of 5% needed to 
explain the replenishment of CRs in our Galaxy. We do not need to assume 
some extreme fraction (several 10%) of energy transfer to CRs in order to 
explain the observed 7-ray flux in particular case. It should be noted that 
by absolute value of energy transfer to galactic CRs, SNRs like Vela Jr. (i.e. 
low energetic SNe) could not power galactic CRs on their own, however the 
lack of energy can be well compensated by the absolute input from more 
energetic than standard 10^^ erg SN explosions. 

We would also like to note that the spectrum in the energy range 0.1-2 
keV is contaminated by the emission of Vela SNR. At the same time the 
emission from Vela Jr. is strongly absorbed in this range. Therefore, it is 
hard to explain the emission of this range that is presented with the light 
grey area at our Figs. [D - Ei 



7. Conclusions 

We applied our hydrodynamical model of the transition stage to Vela Jr. 
Supernova Remnant. We showed that even in moderate ISM density we still 
can account for the considerable 7-ray flux from the SNR. Contrary to the 
widespread opinion that Vela Jr. is a young SNR, our model age is 17 500 yr 
and the remnant is at the beginning of the transition to radiative stage. The 
ISM number density at the place of explosion was 1.5 cm~^ that is normal 
value for the Galactic plane. Explosion energy was 0.2 x 10^^ ergs. The 
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distance to Vela Jr. according to our model is 600 pc. The derived from 
the model broadband spectrum of Vela Jr. as well as corresponding fluxes 
are in good agreement with experimental data. Moreover, the modelled 7- 
ray surface brightness map exhibits features that spatially coincide with the 



nap 

H.E.S.S. 7-ray map in 10(]. The obtained value of the proton cut-off energy 
makes Vela Jr. a probable PeV accelerator. 
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